Wronskian Approach and the One-dimensional Schrodinger 
Equation with Double-well Potential 



Qiu Jian^'^ Ru-Keng Su^'^-'' 
^Department of Physics, Fudan University, Shanghai 200433, P.R.China 
China Center of Advanced Science and Technology (World aboratory), 
P. O. Box 8730, Beijing, P.R.China 
^ qiujian791102@hotmail.com 
Q I ^ rksu@fudan.ac.cn 

o 

(N 

D : Abstract 
O _ 

■ A Wronskian determinant approach is suggested to study the energy and the 
wave function for one-dimensional Schrodinger equation. An integral equation 
and the corresponding Green's function are constructed. As an example, we 
employed this approach to study the problem of double-well potential with 
strong coupling. A series expansion of ground state energy up to the second 

■ order approximation of iterative procedure is given. 
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I. INTRODUCTION 



As a typical bound-state tunnelling problem, the solutions of one-dimensional 
Schrodinger equation with double-well potential are of interest and have attracted much 
attention [1-9]. Recently, Friedberg, Lee and Zhao(FLZ) suggested a new method[7-9] to 
decide the wave functions of Schrodinger equation by quadrature along a single trajectory. 
They gave a good trial wave function for double-well potential and constructed the Green's 
function. By means of the Green's function and the trial wave function, they solved the in- 
tegral equation by iterative procedure and found the successive approximations for the true 
wave function and the corresponding energy. They proved that the iterative solution is con- 
vergent and gave the numerical result of the ground-state energy at first order approximation 
of iterative procedure. 

The objective of this paper is to address another method, namely, the Wronskian deter- 
minant approach(WDA) to study this problem. We will give an integral equation and the 
Green's function by a general argument. In a specific case, our equation reduces to that 
equation given by FLZ method. We will solve our integral equation by iterative procedure 
up to the second order and give the numerical result of the energy. As a comparison, we 
will use the trial wave function given by ref. [9] and the variational method to calculate the 
ground state energy. The numerical results given by the three methods(FLZ, WDA and 
variational method) are quite similar, the differences occur in the higher orders only. This 
result confirms that the Wronskian determinant approach is also successful for studying the 
solutions of one- dimensional Schrodinger equation. 



II. WRONSKIAN DETERMINANT APPROACH 



Consider a second order ordinary differential equation: 





(1) 



with the boundary conditions: 





+ 




(2) 



Assume yi is a solution of the corresponding homogeneous equation: 



y" +Pix)y' + q{x)y = 



(3) 



One can easily prove that another linearly independent solution of Eq. reads: 




(4) 



where A 




is the Wronskian determinant and it satisfies: 




(5) 



2 



Obviously, y2 does not satisfy the boundary conditions, because it must satisfy the unique- 
ness theorem. Suppose the solution y of Eq.(0) can be written as: 

y = C^{x)y^ + C2{x)y2 (6) 
where Ci{x) and C2{x) are two unknown functions which satisfy a constraint: 

C[{x)y^ + C'^{x)y2 = Q (7) 

Since 2/2 does not satisfy the boundary conditions, in order that y satisfy Eq.(||), we have: 

C2(a) = C2{h) = (8) 

Substituting Eq.(|) and Eq.(^ into Eq.(|I|), noting that yi and 7/2 are the solution of Eq.(0), 
we find: 

C[{x)y[ + C',{x)y', = f{x) (9) 
Combining Eq.(0) and Eq.(||), we get: 



Hence: 



yi y2\ (C[\^ (0 
y[y'2 [cU [ fix) 



Ci\ f f yi y2\ ^ f 



C2 J \ y[ yU [fix] 



(10) 



dx (11) 



Ci = - jl\y2f{x)dx (12) 
C2 = l^^yj{x)dx (13) 

y = -yi -^y2f{,x)dx + y2 —yif{x)dx (14) 

where A and B are two constants. From the identities Eq.(^, we have: 

1. The integraton bound A=a 

2. The integration: 

Eq.(|l5D is called the 'resonant condition'. 



The integral bound 'B' is still undetermined. We will discuss this problem in the next 
section. 



III. SCHRODINGER EQUATION 



The one-dimensional Schrodinger equation reads: 

1 

— $ + V<^ = E<^ (16) 

2 ctx^ 

The boundary conditions for bound-states are: 

$(oo) = $(-oo) = (17) 
Consider a perturbation equation: 

1 

--—^ + {V + uj)^ = {E + e)^ (18) 
2 ax^ 

where u is the perturbation potential and e the energy shift. Now we study how to find 
the bound-state wave function and the energy E + e provided $ and E are known. The 
boundary conditions for \I' reads: 

*(oo) = *(-oo) = (19) 

Rewrite Eq.(|I^) as: 



■^-^^ + {V - E)^ = {e - uj)^ (20) 



and take the right hand side of Eq.(^) as the inhomogeneous term. By using Eq.(|l^, we 
have: 

^ = -2$ r<l>i(cu-e)» + 2<l>i r ^uJ-e)^dx (21) 

J B J —CO 

where $i = $ Jq ^~'^dx, according to Eq.(|l) The resonant condition now becomes: 

/oo 
^UJ - e)-^dx = (22) 
-oo 

We can get the energy shift by the identity Eq. (p2D if \E' is obtained. However, the lower 
bound B of the integral in Eq.(pTD is undetermined. Notice that the different choices of B 
will only lead to an extra constant times $ in Eq . (^If) . Then we can fix B=0 and rewrite 
Eq.(pl]) as: 



^/ = c^-2^ ^i{uj - e)^dx + 2^1 <l>(cj - e)*dx (23) 

Jo J~oo 

Combine the two integrations (using Eq. (P^ ) : 

px poo 

^ = c<l>-2<l> ^-^{y)dy ^z){uj - e)-^ {z)dz (24) 

Jo Jy 



There are two ways to determine the constant c: 
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(i) Assume that the perturbation potential is zero at infinity, then = l,and according 
to Eq.(|2|): 

1 = Irl = c - 2 r <^''iy)dy r mico - e)^iz)dz (25) 

<P(OOj JO Jy 

Substitute Eg. (pS]) back into Eg. (p^D , we have: 

^ = $-2$/ <^-^{y)dy <^{z){uj - e)^{z)dz (26) 

J oo Jy 

This is the resuh given by FLZ method[9]. 

(ii) The integral eguation can be solved by iterative method, using \E'o = $ as the zeroth 
approximation, the corresponding energy shift is: 

according to Eg. (p2D . The first approximation of wave function \l'i is given by: 

P X P OO 

^^ = Co$-2<l>/ (!>-^{y)dy ^z){uJ - ei)^{z)dz (28) 

Jo Jy 

Put Co = 1 + A^, where N is chosen to be: 

^ _ 2 JZo <l>'dx <^>-\y)dy <l>{z){u - e,)<l>{z)dz 

we choose cq to satisfy: 

= ($1$) (30) 
which means the modification (\E'i — $) is orthorgonal to $. Similarly, 62 is given by: 

JZo ^^'^lUjdx 

and the corresponding constant Ciis chosen to satisfy: 

(*2|$) = (^l^*) (32) 

and so on. This method of choosing c features that c is different in each order of iterative 
procedure and the modification of each order is orthorgonal to 
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IV. DOUBLE- WELL POTENTIAL 



As an example, we will employ the wronskian determinant approach to discuss the one- 
dimensional Schrodinger equation with double-well potential. Schr odinger equation reads: 

4^^ + ^^V-l)^^ = i^^ (33) 

where g is a. large parameter in the strong coupling case. 

Since the potential is symmetric, we can solve it in the half space x ^ only. So the 
boundary condition for the ground state becomes: 

^(oo) = ^'(0) = (34) 

The trial wave function given by ref.[9] is: 

$ = exp(-|(x - I'fix + 2)) (35) 



The derivative of $ is not zero at the origin. We must do some modification. Let: 

$ + exp(-f ) exp(f (a; - if^x + 2)) for < x <1 

g+l 3 



I $ + ^exp(-^)$ forx>l ^ ^ 



The difference between <l> and $ is of the order exp(— ^).For strong coupling, g ^ 1, 
the difference is exponentially small and can be neglected. We can still choose $ as our trial 
function, the wave function $ satisfies: 

Treat this equation as the unperturbed one, and — (jq^as the perturbative potential cu. 
Employing the Wronsky determinant approach, the integral equation reads: 

pX C OO 

^ = c<l>-2<l>/ ^-^{y)dy ^{z){uj - e)'^{z)dz (38) 

Jo Jy 

The energy shift is given by: 

~ C ^"^dx ^^^^ 
We solve the integral equation Eq.(|38|) by iterative procedure. 

A. First order approximation 

At the first order approximation, \l/o = $,the energy shift ei and the \E'i are: 

^ /o°° -^exp{-'f{x - Ifix + 2))dx 

r (^exp(-f(x-l)2(x + 2))da; ^ ^ 
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$1 = (1 + A^)$ - 2$ / (!>-\y)dy ^z){uJ - e)^{z)dz (41) 

Jl Jy 

where we have put the integration lower bound to be 1 instead of 0. The difference can 
be cancelled by the modification of N. This alteration will bring some convenience in the 
calculation. N is also chosen to satisfy: 

(^^1$) = ($1$) (42) 
For strong coupling, g ^ 1, using the integral formula: 



we obtain a series of 1/g for energy: 

ei = - ^g-' - 0.1660(7-' - 0.2755(7-^.. 

and 

E = 9-\- ^9'' - 0.1660(7-2 - 0.2755^-^.. 
by a straight calculation. 

B. Second order approximation 

At the second order approximation, the energy shift 62 is: 

(1 + N) /o°° ^Mx - 2 /o°° ^Mx /f $-2(y)rfy /,°° ^z){lu - e,)^z)dz 
= l^x ^''^ 



We are facing on the triple integration in the numerator of Eg. (|4^) . The detailed calculation 
of the triple integration is presented in the appendix. The result is: 

1 9 

62 = - 0.1738^-2 - 0.3107^-^.. 

4 64^ y y 



and 



E = g-^- ^g-' - 0.1738(7"' - OMOYg-'.. 
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C. Variational Method 



In this section, we use another method, namely, variational method, to calculate the 
ground state energy. The trial function is chosen as: 

= exp(-^(x - lf{x + 2)) (45) 

where / is the variational parameter. The ground state energy is: 

rtg^F - (TW + W - P){x - 1)^] exp(-f (x - + 2))d^ 

The integration in both the numerator and denominator of Eg. (^Bj) can be calculated by 
using Eg. p3|) . The result is: 

f = g- —g-^ - 0.3320^-2 - 0.8366^^1. . 
64 

and 

1 9 

E = g g-^ - 0.1660^"^ - 0.2855^"^ - 0.6184^"^.. 

4 64 



V. DISCUSSION AND CONCLUSION 

In summary, we have studied the ground state energy and the wave function by using the 
WDA. We have found that the Green's function and the integral eguation are the same as 
that of FLZ method provided the integral lower bound B and parameter c are chosen as oo 
and 1 respectively, we suggest a method to determine parameter c which warrant that the 
modification of $ is orthogonal to $, which is similar to the general perturbation method. 
We have expanded the energy as a series of g~^ for double-well potential and obtained: 

1 9 

E = g — - — '^9^^ " 0.1660(?^^ — 0.2755(7^'^... (1st order iterative approximation) 
1 9 

E = g — - — 0^9^^ ~ 0.17385'^^ — 0.3107^'^^... {2nd order iterative approximation) 
1 9 

E = g — - — ^fl'"^ ~ 0.1660(7~^ — 0.2855(7"^... (variational method) 

The three results are guite similar to each other up to the order 0{g'^) . It is confirmed 
that the WDA is a successful method to investigate the one-dimensional Schrodinger egua- 
tion and the alternative way of determining the parameter c is reasonable. The differences of 
energy between different methods and approximations occur at the order of g~^ and higher 
terms. 
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VI. APPENDIX 

We will calculate the triple integration: 

/ ^uj^dx <^-\y)dy <!>{z){uj ~ ei)<l>{z)dz (47) 

Jo Jl Jy 

Let u := X — 1, the integration becomes: 

oo ru roc 

^uj<^{u)du I ^^^{v)dv ^{w){iu - ei)^{w)dw 

-1 Jo Jv 

/oo rU 
-— — ^ exp{-2g{u^ + u^/?,))du / (2 + vf eM'^div^ + v^/^))dv 
-1 (/ + Uf Jo 

,-^(t^-ei)exp(-25((w^ + wV3))c?w 
V (2 + wY 

oo /-oo /-O /•« /-oo 

...+ /_.../ ... / ... (48) 



10 Jo Jv J-l Jo Jv 

which is the Lagrange's inverse function theorem. 

Suppose function x = f\z) is analytic at the region: |z — zqI < p-T^ + 5 < aig{z — zq) < 
n — 5 : 

f{z) = f{zo) + ak{z - zof + ak+iiz - Zo)''^^ + ... 
while (p{z) is analytic in the neighbourhood of Zo-Then: 

n 1 fjn~l 7 — 7 

^(.) = + E -_|,'(.)(^^^-^-^^)..||„,.(. - (49) 

where Xq = /(^o)- 

When ^{z) = z, Eq. (^9|) gives the formulae of inverse function z = f^^{x). In our case, 
Let: 



X = + u^/3 



We find: 



1/4 -1x^/2 + 11a;... (50) 



(m + 2)2 ' 4 48' 



^ l/16--x'/^ + —X... (51) 



(m + 2)4 ' 8 96 

du={lx-^/^ -l/Q + -^x^/\..)dx (52) 
2 48 

Eqs.([50|)~(|52D are only valid for u > 0. When m < 0, they can be calculated similarly. Thus 
the first part of Eq. (^Sf ) becomes: 



roc 1,1, , 1 , 

/ (1/16- -x'/')i-x-'/' -1/6) expi-2gx)dx (4 + -1/2 _ ^/g) 

JO o 2 Jo 2 

POO ]^ 11 

exp(2^l/)rfy (1/4 - -^V2)(_i/4 + _^i/2 _ ei)(-z-i/2 - 1/6) exp(-2(7;2)rfz (53) 
The problem now is reduced to the calculating of the standard form: 

/ x^'^ exv{-'^gx)dx / y"/'^ exp{2gy)dy / z^^"^ exp{-2gz)dz (54) 

Jo Jo Jy 



where m, n, p are integers ranging from -1 to infinity. 
Rewrite Eq.(p3) as: 



p 00 px p 00 

/ x""^^ exp{-2gx)dx y""^^ exp{2gy)dy z^'^ exp{-2gz)dz + 

Jo Jo Jx 

P CXD P X PX 

/ x""^^ exp{-2gx)dx / y""'^ exp{2gy)dy / z^^^ exp{-2gz)dz (55) 

Jo Jo Jv 



In the first part, let: x = zs,y = zst, z = z s G [0, 1] t G [0, 1] and J = 

In the second part, let: x = x,y = xst, z = xs s G [0, 1] t G [0, 1] and J 
Eq.( becomes: 



d(x,y,z) 
d{z,s,t) 
x,y,z) 



d(x,s,t) 



Z^S 
2 

= X S 



where 



/ s-^f'l^z s exp \-2gz (1 - + s)\ dzdsdt + (56) 

s 2 t / X 2 exp [-25(x [1 - st + s)] dxdsdt 



: 2 G [0,cx)),s G [0, l],t G [0,1] (57) 
: a; G [0,cx)),s G [0, l],t G [0, 1] (58) 



The variables x and z can be integrated Eg . (|56|) becomes: 

r/ m+n+p+6 \ 
1 r5 I / m + n+2 ,r, m+n+p+4 m + n+p+6 

, J.,^ / s^^r/^z — 2 1 _ St + s 2 c^srft + 59) 

-f^/ m+n+p+6 \ 

i I ;^ I / m+p+2 ,„ m + n+p+4 m + n+p+6 

s^^r/\ — 5 — (i- st + s) 2 (^s(^t 



(2^)(™+"+P+6)/2 

Now the parameter g is extracted out, the remaining two double integrations are inde- 
pendent of g and the area for integration is a rectangle, they can be calculated easily. 

We calculate the integration which satisfy ^ g^^id sum all these integrations, 

it will give us the first part of Eq.(^). The second part can be calculated with the same 
procedure. 

We thank Prof. W.Q.Zhao and Prof. S.Z.Hu for helpful discussions, this work was 
supported in part by NNSF of China under contracts No.l9975050, 10047045, 19947001. 
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